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Jtethod for designing the spider spring of a control cluster 
of a nuclear fuel assem bly, a corresponding system, computer 

programme and product 

The present invention relates to a method for designing a 
nuclear fuel assembly which is intended to be positioned in 
a nuclear reactor, the assembly comprising a plurality of 
guide tubes,- and a control cluster which itself comprises a 
plurality of control rods which are received in the guide 
tubes and a support for control rods, the assembly compris- 
ing a helical spring for damping the impact of the support 
against an upper end piece of the assembly in the event of 
the control cluster falling during a shutdown of the nuclear 
reactor . 

It will, be appreciated that nuclear fuel assemblies must be 
dependable in order to allow reliable operation of nuclear 
reactors i 

* 

Thus, design and construction provisions for such assemblies 
have been drawn up. 

These provisions impose a general framework and minimum cri- 
teria which the assembly constructors must take into consid- 
eration. 

As far as the helical damping spring is concerned, the de- 
sign provisions require verification by means of tests that 
the integrity of the spring has not been affected during the 
impact brought about in the event of a shutdown of the reac- 
tor . 

Although the criterion imposed by the design provisions al- 
lows assemblies to be designed with satisfactory reliabil- 
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ity, it would be desirable to limit the safety margins dur- 
ing design in order to reduce the mass and the cost of the 
assemblies constructed. 

An object of the invention is to overcome this problem by 
providing a method which allows reliable nuclear fuel assem- 
blies to be designed, whilst limiting the design margins. 

To this end, the invention relates to a method for designing 
a nuclear fuel assembly which is intended to be positioned 
in a nuclear reactor, the assembly comprising a plurality of 
guide tubes and a control cluster which itself comprises a 
plurality of control rods which are received in the guide 
tubes and a support for control rods, the assembly compris- 
ing a helical spring for damping the impact of the support 
against an upper end piece of the assembly in the event of 
the control cluster falling during a shutdown of the nuclear 
reactor, characterised in that the method comprises, the 
following steps : 

a) establishing the progression of the speed of the con- 
trol cluster after the impact of the support against the up- 
per end piece, 

b) establishing, based on the speed established in step 
a) , a maximum longitudinal load for compression of the 
spring, and 

c) establishing, based on the maximum longitudinal load 
for compression, at least a maximum shearing stress in the 
spring . 

According to specific embodiments, the method can comprise 
one or more of the following features, taken in isolation or 
according to all technically feasible combinations: 
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- a. maximum shearing stress is a shearing stress along the 
neutral axis of the spring, 

- a maximum shearing stress is a shearing stress along the 
axis of the spring nearest the longitudinal centre axis 
thereof , 

- the method further comprises a step for verifying, using a 
maximum shearing stress established in step c) , that a maxi- 
mum stress admissible by the spring has not been exceeded. 

The invention further relates to a system for designing a 
nuclear fuel assembly, characterised in that it comprises 
means for carrying out the steps of a method as defined 
above . 

According to a variant, the system comprises a computer and 
storage means, in which at least a programme comprising in- 
structions for carrying out steps of the method for design- 
ing a nuclear fuel assembly is stored. 

The invention further relates to a computer programme com- 
prising instructions for carrying out the steps of a method 
as defined above. 

The invention also relates to a medium which can be used in 
a computer and on which a programme as defined above is re- 
corded . 

The invention will be better understood from a reading of 
the description below which is given purely by way of exam- 
ple with reference to the appended drawings, in which: 

- Figure 1 is a schematic, perspective cut-away view of a 
nuclear fuel assembly which is designed by a method accord- 
ing to the invention, 
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- Figure 2 is a schematic, partially sectioned side view 
drawn to an enlarged scale of the structure of the spider of 
the assembly of Figure 1, 

- Figure 3 is a partial schematic side view of the assembly 
of Figure 1, illustrating more particularly a pair compris- 
ing a guide tube/control rod, 

- Figure 4 is a block diagram illustrating the system for 
designing the assembly of Figure 1, 

- Figure 5 is a flow chart illustrating successive steps of 
the design method carried out by the system of Figure 4, 

- Figure 6 is a progression curve of the falling speed of a 
control rod before it is introduced in the lower portion of 
the corresponding guide tube, this progression being calcu- 
lated by the system of Figure 4, and 

- Figure 7 is a progression curve of the falling speed of 
the same control rod in the lower portion of the correspond- 
ing guide tube, this progression being calculated by the 
system of Figure 4. 

Figure 1 illustrates a nuclear fuel assembly l which mainly 
comprises a square-based lattice 2 for nuclear fuel rods 3 
and a control cluster 4 . 

The assembly 1 comprises grids 5 for maintaining the rods 3, 
which grids 5 are distributed over the height of the rods 3 . 
A lower end piece 6 is arranged under the lower ends of the 
rods 3 and an upper end piece 7 above the upper ends of the 
.rods 3. The upper end piece 7 is provided with springs 8 for 
pressing against the upper bearing plate of the reactor 
core, in which the assembly 1 is intended to be placed. 

The control cluster 4 comprises a plurality of control rods 
10, for example, 24. Conventionally, the control rods 10 
comprise a material which absorbs neutrons. 
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The rods 3 and 10 extend . in parallel with a vertical longi- 
tudinal direction L. 

The rods 10 are carried at the upper ends thereof by a sup- 
port ll which is generally referred to as a spider. 

As illustrated more particularly in Figure 2, the spider 11 
comprises a vertical central upper head 12 and a series of 
arms or vanes 13 which extend radially outwards from the 
lower end of the upper head 12 as far as the radially outer 
ends 14 thereof. 

Each control rod 10 is connected to an arm 13 at the upper 
end thereof. 

The upper head 12 of the spider 11 has a central blind hole 
15 which opens towards the bottom and in which a damping 
helical spring 16 is received. The spring 16 extends verti- 
cally along a centre axis A. A tightening screw 17 extends 
substantially over the entire height of the hole 15 and is 
screwed into the wall 18 delimiting the upper portion of the 
hole 15. 

The lower portion of the screw 17 extends through the base 
of. a retaining ring 20 which rests on the lower end of the 
spring 16. The head 21 of the screw 17 rests, at the top, 
against the base of the retaining ring 20 in order to press 
the spring 16 against the wall 18 of the upper head 12. 

As illustrated in Figure 3 for a control rod 10, each con- 
trol rod 10 is received in a respective guide tube 24 which 
is arranged in the lattice 2 of fuel rods 3. In this manner, 
24 pairs comprising a guide tube/control rod are formed. 
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Since each of these pairs has a similar structure, only one 
will be described below . 

The guide tube 24 extends from the lower end piece 6 a_s far 
as the upper end piece 7. The guide tube 24 comprises a 
lower portion 26 of reduced inside diameter and an upper 
portion 27. The lower portion 26 is connected to the lower 
end piece 6 by a collared screw 28, through which a vertical 
through-hole 29 extends. 

The lower portion 26 of the guide tube 24 surrounds the con- 
trol rod 10 with a radial passage gap J. 

The upper portion 27 is fixed to the upper end piece 7 and 
opens at the outside of the assembly 1, 

Lateral apertures 30, only one of which can be seen in Fig- 
ure 4, are provided in the upper portion 27 near the lower 
portion 26. 

When the assembly 1 is placed in a nuclear reactor, the 
cooling liquid of the reactor fills the interior of the 
guide tube 24. 

Conventionally, the control cluster 4 can be moved verti- 
cally relative to the remainder of the assembly 1 in order 
to allow adjustment of the reactivity during normal opera- 
tion of the reactor, and therefore variations in power from 
zero power up to maximum output depending on the vertical 
introduction of the control rods 10 in the lattice 2 of rods 
3 . The vertical displacement of the control cluster 24 is 
conventionally carried out by way of a drive rod which is 
connected to the upper end of the upper head 12 . 
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When the reactor is shut down, the drive rod and the assem- 
bly 4 fall owing to gravity. 

At the start of this falling movement, the control rods 10 
are guided only by the upper portions 27 of the guide tubes 
24 and have not yet reached the lower portions 26. 

Once the falling action has ended, the lower ends of the 
control rods 10 are introduced in the lower portions 26. The 
cooling fluid contained in the portions 26 is then violently 
forced, on the one hand, upwards thereby and, on the other 
hand, downwards through the apertures 29 of the collared 
screws 28. 

Each lower portion 26 therefore behaves in the manner of a 
hydraulic damper braking the falling movement of the corre- 
sponding control rod 10, and therefore of the assembly 4. 

This braking phase ends at the end of the travel path with 
the impact of the spider 11 against the upper end piece 7 of 
the assembly 1 . 

This impact is carried out by means of the retaining ring 
20. During this impact, the spring 16 is compressed verti- 
cally in order to absorb the shock. 

According to the invention, the assembly 1 has been designed 
in order to take into consideration the specific stresses 
brought about in the assembly by the fall of the control 
cluster 4 during such a shutdown of the reactor. 

In this manner, in order to design the assembly 1, in par- 
ticular a data-processing system 32 has been used, as illus- 
trated schematically in Figure 4. 
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This system 32 comprises, for example, a computer or data 
processing unit 34 comprising one or more processors, stor- 
age means 36 , input/output means 38, and optionally display, 
means 4 0 . 

Instructions which can be carried out by the computer 34 are 
stored in the form of one or more programs in the storage 
means 36. 

These instructions are, for example, instructions in FORTRAN 
programming code. 

These various instructions, when they are carried out by the 
computer 3.4, allow the method illustrated by the flow chart 
of Figure 5 to be carried out. 

In a first step illustrated by the box 42 of this Figure, 
the computer 34 calculates, based on data 43, the develop- 
ment of the falling speed of a control rod 10 in the upper 
portion 27 of the corresponding guide tube 24 in the event 
of a shutdown of the reactor. 

This calculation can be carried out assuming, for example, 
that the control rod 10 is first subjected to constant 
loads : 

gravitational force: fg « Mg, 

Archimedes' thrust: fa = -p gv, 

pressure difference in the core: fc, and 

mechanical friction: fm, 
where M and V are the mass and the volume, respectively, of 
the assembly 4 and the drive rod thereof. 
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The control rod 10 is also subjected to loads as a function 
of the speed or position thereof, for example, hydraulic 
friction which can be obtained from: fh = -cl (M + pV) v 2 , 
with v = speed of the assembly 4 and therefore of the rod 10 
in question. 

Thus, the equation of the movement of the rod in the upper 
portion 27 of the guide tube 24 is as follows; 

(M*pV)^-zf 



This gives: 



with cl = hydraulic friction in the guide tube, and 

c2 - fs+fa+fc+fm 
M+pV 

Cl and c2 are, for example, experimental data measured dur- 
ing drop tests of the control cluster 4. These data are, 
with the other data necessary for the calculation, such as 
the mass and the volume of the assembly 4 and the drive rod 
thereof, introduced, for example, in the form of a file 43 
by way of the input /output means 38. 

The computer 34 resolves the equation of the movement of the 
control rod 10, for example, using the NEWTON method. 

Thus, the progression of the speed of the control rod 10 in 
the upper portion 27 is known as a function of time. The 
profile established in this manner can be displayed in the 
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form of a curve by the display means 40. This curve is il- 
lustrated by Figure 6 . 

In this manner, at the end of the step illustrated by the 
box 42, the speed of the control rod 10 is known at tfcie 
point of entry to the lower damping portion 26 of the guide 
tube 24 . 

Based on the results of the step of box 42, the computer 34 
calculates the progression of the speed of the control rod 
10 during its fall in the lower damping portion 26. 

This step is schematically illustrated by the box 4=4- 

This step can be carried out using the following equation: 

with c2= /gV* = M ~P V „ 

M+pV M+pV g 

SCA = cross -section of the rod 10 and 

NCA *= number of rods 10 in the assembly 4. 

Therefore, the hypothesis that f G and f m are negligible is 
applied here. 

The difference AP represents the elevated pressure produced 
in the cooling liquid contained in the guide tube 24, that 
is to say, the pressure thereof between the lower end of the 
rod 10 and the pressure present in the upper portion 27 of 
the guide tube 24 . 

AP can be established by the following formula: 




y M+JpV 
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Ap = yPQ 2 v 2 (EXPA+CONTRA+FECRjcCISAX2) 

where EXPA = U^jA M CM ]f 
VACM* SACTG " 

with SM = cross-section of the lower portion 26, 

SACM = SM - SCA •» cross -section of the annular space between 

the rod 10 and the lower portion 26, 

SACTG = STG - SCA, where STG is the cross -section of the up- 
per portion 27 of the guide tube 24, 

CONTRA =0.4 (t-^d£^\_SCA_) 
" SM *SACM> ' 

PECR = coefficient of loss of load owing to friction in the 
lower portion 26, 

CISA = 1 

DM = mean diameter of the guide tube 24 in the upper portion 
27, 

z = height of the rod 10 introduced in the lower portion 26 
of the guide tube 24, and 

Q = fraction of liquid flowing upwards out of the lower por- 
tion 26 . 

The resolution of the equations governing the movement of 
the rod 10 after entry into the lower portion 26 is carried 
out by the computer 34, for example, using the RUNGE - KUTTA 
method. 

Thus, at the end of the step 44, the progression of the 
speed of the control rod 10 in the lower portion 26 of the 
guide tube 24 is known before the impact of the spider 11 on 
the upper end piece 7. 

The speed profile established in this manner can be dis- 
played, for example, by the means 40. as illustrated in Fig- 
ure 7. on the curve in Figure 7, the speed profile estab- 
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Lished during step 44 is the portion located to the left of 
the point 45 • 

The computer 34 then carries out, in the step of box 46, the 
calculation of the maximum elevated pressure produced APmax- 

This calculation can be carried out, for example, based on 
the formula : 

AP - i pQ^iEXPA^ONTRA^FECFxClSAxz) . 

The computer 34 carries out, in the step of box 48, -the cal- 
culation of a circumferential stress and maximum normal aewoc# 
to which the lower portion 26 of the guide tube 24 is sub- 
jected owing to the maximum elevated pressure APmax. 

This stress can be calculated based on the formula: 



where DPM = inside diameter of the lower portion 26 and 

EMP » minimum thickness of the wall of the lower portion 26 . 

The system 32 can then provide, owing to the input /output 
means 38, a first result in the form of a file 49 containing 
the value Oomax established, and optionally the maximum ele- 
vated pressure APmax established. 

Next, the system 32 carries out the calculation of the pro- 
gression of the speed of the control rod 10 after it comes 
into contact with the spider 11 and the upper end piece 7 . 




Fax emis par: 33 1 S3 20 14 91 CABINET IAVOIX le 26/08/04 17:37 A4 NORM Pg: 19/28 



13 

This calculation step is illustrated by the box 50 in Figure 

This calculation can be carried out, for example, using the 
following equation when the ring 20, and therefore the spi- 
der 11, is in contact with the upper end piece 7: 

(M + pV) ^ =(M- P V)g-PRCH-K(z-LAJ)-c3v 

with PRCH = pretension of the spring 16 * PRCMP x K, where 
PRCMP is the precompression of the spring 16 and K the ri- 
gidity of the spring 16, 

LAI n distance travelled by the control rod in the lower 
portion 26 before impact, and 

c3 = coefficient of hydraulic damping in order to model the 
damping in the lower portion 26. 

In the event of a rebound, that is to say, when the spider 
11 is no longer in contact with the upper end piece 7, the 
equation for movement of the control rod 10 in question is 
written as follows: 

(M + p v) = (M - P V) g - c3 v 

These two equations are integrated by the computer 34, for 
example, using the RUNGE-KUTTA method. 

Therefore, the step 50 allows the kinematics of the control 
cluster 4 to be established during the mechanical damping of 
the shock by the spring 16. The speed profile established in 
this manner can be displayed, for example, by the means 40. 
This profile corresponds to the portion located to the right 
of the point 4 5 on the curve in Figure 7. 

Based on the results of this step, the system 3 2 carries 
out, in the step 52, the calculation of a maximum vertical 
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compression force Fmax, to which the spring 16 is subjected 
during the mechanical damping. 

This calculation can be carried out, for example, based on 
the following formula: 

Fmax = max {K(z-LAI) + PRCH} 
The system 32 then carries out, in the step of box 54, the 
calculation of an approximate maximum shearing stress tmax i^ 
the spring 16 : 

^ _ SFmaxDFN 

T MAX — f= , 

tuDFR 

with DFN = DER-DFR and 

DER = outside diameter of the spring 16, 
DFR = diameter of the wire of the spring 16. 

Subsequently, the system 32 can optionally carry out, based 
on the maximum stress Tmax, the calculation of maximum cor- 
rected stresses-. 

These stresses can be calculated by multiplying Tmax by dif- 
ferent factors. 

Thus, it is possible to calculate: 

xmaxi = tmax X Kc, and 

with Kc= 1 + 
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The stress Tmaxi corresponds to the shearing stress along the 
neutral axis FN (Figure 2) of the spring 16. The stress Tmax2 
corresponds to the stress along the axis F2 (Figure 2) of 
the spring 16 nearest the vertical centre axis A of the 
spring 16 (see Figure 2) . 

At the end of this step illustrated by the box 56, the sys- 
tem 32 provides the various maximum shearing stresses calcu- 
lated, for example, in the form of data stored in a file 57 , 
which are transmitted by the input/output means 38. 

Based on the data contained in the files 49 and 57, which 
have also been stored in the storage means 36, the "computer 
34 will verify that the maximum stresses calculated are in- 
deed acceptable for the materials which respectively consti- 
tute the guide tube 24 and the helical spring 16- 

This step has been schematically illustrated by the box 58 
in Figure 5. During such a step, the system 32 will, for ex- 
ample, verify that the maximum shearing stresses calculated 
during the steps 54 and 56 are less than maximum values ad- 
missible by the material which constitutes the spring 16. 
This verification is carried out by a comparison of Tmax* Tmaxi 
and tiuoa.with a maximum value admissible by the material of 
the spring 16. 

As far as the maximum circumferential stress Gqkax is con- 
cerned, the verification can be carried out based on a for- 
mula of the type : 

admissible 

where a a en,issibie refers to the material which constitutes the 
lower portions 26 of the guide tubes 24. 
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The function f can be a function which takes into considera- 
tion other stresses to which the guide tubes 24 can be sub- 
jected. Such a stress can be a vertical compression stress 
ct a , to which the guide tubes 24 are subjected during the con- 
tact of the springs 8 of the upper end piece 7 against the 
upper bearing plate of the core in order to counterbalance 
the hydrostatic thrust during operation. 

Thus, the function f can be, for example, in the form of f 

(<J0MAX, CJ A ) = CF0MAX + CT A . 

It will be appreciated that this last step, illustrated by 
the box 58, can be carried out by separate software which 
generally carries out the validation of various design pa- 
rameters of the assembly 1 based on results provided by 
various pieces of software each dedicated to taking into 
consideration specific operating conditions and which in- 
clude the software which carries out the steps 42, 44, 46, 
48, 50, 52, 54 and 56. 

In general terms, the file 43 comprising the data 43 used by 
the method for the various calculations can comprise the 
data of Table l below. 
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outside diameter of control rod 10 

inside diameter of upper portion 27 

inside diameter of lower portion 26 

total length of lower portion 26 

damping travel before impact 

minimum thickness of wall of lower portion 26 

maximum roughness of rod 10/tube 24 

diameter of aperture 29 

length of aperture 29 

roughness of aperture 29 

moving mass M 

volumetric mass of liquid 

kinematic viscosity of liquid 

cl 

c2 

Young's modulus of guide tube 24 

Poisson's ratio of guide tube 24 

spring precompression 16 

preloading of spring 1 6 

length of spring 16 with contiguous turns 

outside diameter of spring 1 6 

diameter of wire of spring 16 

compression when upper head 12 is in contact with 

upper end piece 7. 

END _ 



(m) 
(m) 
(m) 
(m) 
(m) 
(m) 
(m) 
(m) 
(m) 
(m) 

(k^m 3 ) 

(n?/s) 

(/m) 

(m/s 2 ) 

(Pa) 

(m) 
(N) 
(m) 
(m) 
(m) 
(m) 



Nominal; maximum 
Nominal; maximum 
Nominal; maximum 



Table 1 



Similarly, the file 49 con^prising the results from step 48 
can comprise the data of Table 2 below • 



APmax: maximum elevated pressure in lower portion 26 

Zmax* corresponding penetration in lower portion 26 

ogmax" maximum stress in lower portion 26 

ftnax: maximum force on lower end piece 6 

tdun duration of fall in lower portion 26 before impact 

vfui: speed of impact of assembly 4 on upper end piece 7 



(Pa) 
(m) 
(Pa) 
(N) 

00 
(m/s) 



Table 2 



The file 57 comprising the results of step 56 can itself 
contain the data of Table 3 below. 



Fmax : maximum compression force on spring 16 
hivtAx: maximum deflection of spring 1 6 



<N) 
(in) 
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xmax: approximate mavimum stress in spring 
^maxi * approximate maximum stress corrected by Kc 
*max2: approximate maximum stress corrected by K (Wahl coef- 
ficient) 



(Pa) 
(Pa) 
(Pa) 



Table 3 



It has been possible to verify by experiment that the maxi- 
mum elevated pressures and the maximum stresses obtained by 
means of steps 42 , 44, 46 and 48 were reliable. In this man- 
ner, the first corresponding part of the method allows reli- 
able guide tubes 24 to be designed. Furthermore, this first 
part calculates only a single stress which appears to be the 
pertinent stress for the conditions being considered* Conse- 
quently, this first part of the method allows the security 
margins to be limited during design, and therefore assem- 
blies which are relatively light and economical to be de- 
signed. 

The second part of the method, which corresponds to steps 
50,-52, 54 and 56, also allows maximum stresses to be relia- 
bly calculated, as confirmed by experiment. 

Thus, the second part of the method allows a reliable design 
to be arrived at by calculation for the spider springs 16, 
which design is found to be advantageous in comparison with 
the method of tests alone which is currently imposed by pro- 
visions. It will be appreciated that the second part of the 
method calculates only the small number of stresses, and in 
particular those on the axis F2 of the spring 16 nearest the 
centre axis A of the spring, which are found to be pertinent 
to the conditions envisaged. In this manner, the second part 
of the method allows the design margins to be reduced . 
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In more general terms, the steps 42, 44 , 46 and 48, on the 
one hand, and 50, 52, 54 and 56, on the other, can be car- 
ried out by separate pieces of software. 

In order to increase the reliability of the calculation, for 
carrying out the first part of the method it is possible to 
use, as the passage gap J, the nominal value of the gap, or 
this nominal value corrected by the manufacturing tolerance, 
or a value resulting from statistical studies of the distri- 
bution of passage gaps J obtained in constructed assemblies. 

In a variant, it is possible to use a gap value J which is 
greater for steps 42 and 44 and a smaller gap value J for 
steps 46 and 48. This allows a high stress value aeMwc to be 
calculated because the speed reached during the fall of the 
rod 10 in question is high and the volume available in the 
lower portion 26 for the liquid during damping is small. 
However, this high stress value is not unrealistic and 
therefore does not lead to unjustified design margins, as 
illustrated by the following example. 

According to a specific variant, the upper value can be a 
maximum value for gap J which is verified with a given prob- 
ability, for example, 95%, in constructed assemblies, and 
the lower value can be a minimum value obtained with the 
same probability. This variant allows an approximation of 
the situation where a single pair comprising a guide 
tube/control rod has minimum gap J, where the maximum stress 
Oomax would be reached, and where all the other pairs compris- 
ing a guide tube/control rod have the maximum passage gap or, 
which would be the most extreme case. 

In some variants, the first part of the method could also 
take into consideration forms of the lower damping portion 
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26 which are different from those described previously. In 
this manner, these lower damping portions could have a plu- 
rality of successive portions of reduced diameter, option- 
ally separated by portions of increased diameter, generally 
referred to as cavities. In some variants, the first part of 
the method is carried out with collared screws 28 which are 
not perforated by holes 29. 

In still more general terms, the first part and the second 
part of the design method described can be used independ- 
ently of each othe:r. In this manner, it is possible to carry 
out the second part relating to the design of the spring 16 
without referring to the calculation of the elevated pres- 
sure AP and the stress crew*. 
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